Prawdopodobnie rozkład długości hospitalizacji nie zależy od charakteru udaru, co oznacza, że powtórny udar leczy się w ten sam sposób (stosując te same procedury)
zgony <- paste0("pZgon", c(7, 30, 90, 365))
wiek_plec_m <- dane[!(transport) & PLEC == "M" & lUdar < 3 & lUdar == kUdar,
.(pZgon7 = 100*mean(zgon7),
pZgon30 = 100*mean(zgon30),
pZgon90 = 100*mean(zgon90),
pZgon365 = 100*mean(zgon365)
),
by = c("WIEK", "kUdar")
]
wiek_plec_k <- dane[!(transport) & PLEC == "K" & lUdar < 3 & lUdar == kUdar,
.(pZgon7 = 100*mean(zgon7),
pZgon30 = 100*mean(zgon30),
pZgon90 = 100*mean(zgon90),
pZgon365 = 100*mean(zgon365)
),
by = c("WIEK", "kUdar")]
# plot_ly(wiek_plec_7_m,
# x=~WIEK, y=~pZgon7,
# color = ~kUdar,
# colors = kolory,
# type = "bar") %>% layout(barmode = "overlay",
# bargap = .1,
# hovermode = "x"
# )
#
# plot_ly(wiek_plec_7_k,
# x=~WIEK,
# y=~pZgon7,
# color = ~kUdar,
# colors = kolory,
# type = "bar") %>% layout(barmode = "overlay",
# bargap = .1,
# hovermode = "x"
# )
subplot(lapply(zgony, function(x) plot_ly(wiek_plec_k,
x=~WIEK,
y=~get(x, wiek_plec_k),
color = ~kUdar,
colors = kolory,
showlegend = "if"(x == zgony[1], FALSE),
type = "bar") %>% layout(annotations = list(x = 0.2 ,
y = 1.05,
text = paste0(letters[which(zgony == x)],
") śmiertelność ",
stri_extract(x,
regex = "[[:digit:]]*$"),
"-dniowa") ,
showarrow = FALSE,
xref="paper",
yref="paper"
),
title = "Śmiertelność kobiet w zależności od wieku",
barmode = "overlay",
bargap = .1,
hovermode = "x")
), nrows = 2)
subplot(lapply(zgony, function(x) plot_ly(wiek_plec_m,
x=~WIEK,
y=~get(x, wiek_plec_m),
color = ~kUdar,
colors = kolory,
showlegend = "if"(x == zgony[1], FALSE),
type = "bar") %>% layout(annotations = list(x = 0.2 ,
y = 1.05,
text = paste0(letters[which(zgony == x)],
") śmiertelność ",
stri_extract(x,
regex = "[[:digit:]]*$"),
"-dniowa") ,
showarrow = FALSE,
xref="paper",
yref="paper"
),
title = "Śmiertelność kobiet w zależności od wieku",
barmode = "overlay",
bargap = .1,
hovermode = "x")
), nrows = 2, margin = .1)
Powtórny udar dotyczy węższej grupy wiekowej. W grupie najstarszych pacjentów drugi udar nie zdarza się ze względu na wysoką śmiertelność. Można spróbować też dla grup wiekowych (2 różne podziały)
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in `[.data.table`(popWiekBK_M, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 1 column 'lPacj': 18.961507 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in `[.data.table`(popWiekBK_M, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 2 column 'lPacj': 12.561175 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in `[.data.table`(popWiekBK_K, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 1 column 'lPacj': 18.961507 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in `[.data.table`(popWiekBK_K, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 2 column 'lPacj': 12.561175 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in `[.data.table`(popWiekEks_M, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 1 column 'lPacj': 33.072466 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in `[.data.table`(popWiekEks_M, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 2 column 'lPacj': 4.642784 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in `[.data.table`(popWiekEks_K, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 1 column 'lPacj': 33.072466 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in `[.data.table`(popWiekEks_K, , `:=`(lPacj, 100 * lPacj/
## sum(lPacj)), : Group 2 column 'lPacj': 4.642784 (type 'double') at RHS
## position 1 truncated (precision lost) when assigning to type 'integer'
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
## Warning in RColorBrewer::brewer.pal(N, "Set2"): minimal value for n is 3, returning requested palette with 3 different levels
To jest do omówienia z ekspertami i z prof. Kamińskim.